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Decoherence-Free Subsystems (DPS) are a powerful means of protecting quantum information 
against noise with known symmetry properties. Although Hamiltonians theoretically exist that can 
implement a universal set of logic gates on DFS encoded qubits without ever leaving the protected 
subsystem, the natural Hamiltonians that are available in specific implementations do not necessarily 
have this property. Here we describe some of the principles that can be used in such cases to 
operate on encoded qubits without losing the protection offered by the DFS. In particular, we show 
how dynamical decoupling can be used to control decoherence during the unavoidable excursions 
outside of the DFS. By means of cumulant expansions, we show how the fidelity of quantum gates 
ja. ' implemented by this method on a simple two-physical-qubit DFS depends on the correlation time 

^--s \ of the noise responsible for decoherence. We further show by means of numerical simulations how 

^vj . our previously introduced "strongly modulating pulses" for NMR quantum information processing 

' can permit high-fidelity operations on multiple DFS encoded qubits in practice, provided that the 

f) ' rate at which the system can be modulated is fast compared to the correlation time of the noise. 

^j^ , The principles thereby illustrated are expected to be broadly applicable to many implementations 

^^ ■ of quantum information processors based on DFS encoded qubits. 

(N. 

I. INTRODUCTION 

> ■ . . 

*ys ' The interaction of a quantum system with its environment leads to the loss of quantum phase information and 

f~*) . interference, which is widely known as decoherence [llQ,y{. This conversion of quantum information into classical 
^NJ ' information is a major obstacle in any application where pure quantum phenomena are sought, particularly quantum 
^^ , information processing (QIP). One promising strategy for countering decoherence is to encode the information in 
^^ ■ subsystems that, because of the symmetry of the interactions between the qubits and their environment, arc invariant 
^^ ' under the action of the noise generators H, [3 . 

^r • Universal computation within these decoherence-free subsystems (DFS) has been shown to be possible from a 

1-^ '■ theoretical point of view |a, Qilal- ^^ particular, universal fault-tolerant computation within a DFS is possible if the 
^H. exchange interaction between qubits can be switched off and on at will [3, [ifl HH ■ The issue that we address here 
-j-J ' arises in systems where this Hamiltonian does not occur naturally. 

^ , The total Hamiltonian is conveniently divided into a time-independent part. Hint, and a part that depends on a set of 

1-^ ' experimentally controllable, time-dependent parameters, Hextiptii^)) ■ There are many useful DFS encodings for which 
qh] the generators of logical qubit rotations are not contained in the set of available total Hamiltonians: Htot = Hint + Hext- 
• • ■ As long as the system is universal, any desired propagator on the logical subsystem can be composed from the evolution 
.5^ ' under a series of time dependent external Hamiltonians; the instantaneous total Hamiltonian need not preserve the 
S^ i protected subsystem. In such cases, the extent to which universal fault-tolerant computation is possible depends on 
^ ' the details of the control fields, as well as the spectral density of the noise ^^. This paper will discuss these issues in 
the context of a specific QIP implementation based on liquid-state NMR. 

In Section II it is shown that the radio-frequency (RF) control fields used in NMR necessarily cause the encoded 
information to "leak" from the protected subsystem into other parts of the total system Hilbcrt space, where it is 
subject to decoherence. This is illustrated by numerical simulation of a simple example - the encoding of one logical 
qubit in a two-spin decoherence-free subspace for collective dephasing |l3l | - which is referred to throughout the paper. 
It is further shown that in the case of two DFS qubits, each with the same noise model, the leakage rate is generally 
nonzero even in the absence of the control Hamiltonian. 

In Section III, we briefly review the application of stochastic Liouville theory [ij, [13 to understand the effective 
decoherence in the presence of an external Hamiltonian that modulates the spin dynamics, as is the case for strongly- 
modulating pulses (SMP) [la. [l7l| or optimal control theory [ig, dynamical decoupling [13 or bang-bang control 
|20l | . We then use these results in Section IV to quantitatively understand the Carr-Purcell [21| (CP) sequence. 
CP is perhaps the original dynamical decoupling sequence and the archtypc for observing the influence of the noise 
correlation time on the effective decoherence rate. 
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We will show that the CP sequence can be effective at suppressing both leakage and decoherence provided that 
one can modulate the system on a time scale shorter than the correlation time of the noise. Under the assumption 
of instantaneous and perfectly selective single spin 7r-pulses, the exact dependence of the overall gate fidelity on the 
pulse rate and the correlation time of the noise is derived for a single DFS qubit. Finally in Section V, we apply these 
ideas to a physical system with realizable models of control fields. This includes limitations on the available RF power 
and the lack of frequency selectivity among the physical spins. As expected, pulses of finite duration degrade the gate 
fidelities of these operations. Analytical solutions are generally not feasible, and hence the amount of degradation 
that can be expected for a range of experimentally realistic parameters is evaluated by means of numerical simulations 
for simple quantum gates operating upon one or two encoded qubits. Lastly, we show the criteria of noise correlation 
times that are compatible with high fidelity control. 

II. LEAKAGE FROM A TWO-LOGICAL QUBIT DFS 

In physical implementations of QIP, qubits arc not embodied in well-isolated two-state systems: rather, they are 
embedded in larger Hilbert spaces containing additional states that arc intended to not participate in the computation. 
In addition, when physical two-state systems are combined into logical qubits for noise protection and/or correction, 
additional redundant degrees of freedom are introduced. Leakage to these "external" degrees of freedom can destroy 
the coherent dynamics of the qubit [22, 123 . 

Such leakage can be introduced by the control fields applied to implement specific logic gates, and even in absence 
of external fields, by the internal Hamiltonian itself. Many modulation methods have been engineered to refocus the 
terms in the internal Hamiltonian responsible for leakage [2J| . We explore the fidelity that can be reasonable expected 
based on the details of the modulation scheme and the spectral density of the noise. 

In this section we use a simple well-known DFS to motivate our discussion. The DFS encodes one logical qubit in 
two physical spin-^ particles and protects against collective dephasing caused by fully correlated uni-axial noise. In 
NMR, for example, fluctuations of the quantizing magnetic field Bz at a local molecule appear fully correlated, yet 
lead to dephasing when averaged over the spin ensemble |l3 | . These fiuctuations are described by a Hamiltonian of 
the form HsE{t) ~ lBz{t) Z , where Z = ^ ^^ (j\ is the total angular momentum of the spins along the z-axis and 
7 is their gyromagnctic ratio. The DFS is based on the encoding |0)l = |01}, \1)l = |10}. A basis for the space of 
operators on the encoded qubit, in turn, is given by the four logical Pauli operators: 

a^ ^ ^{al-a^) o^ ^ \{olal + oy^ (1) 

This two spin-i particle Hilbert space (C"' = C^ ® C^) can be described as a direct-sum of the total angular 
momentum subspaces, Zq © Zj^i © 2^_i, where I is the total angular momentum projected along the quantization axis. 
The logical basis states |0)l and |1)l reside exclusively in Zq, where Zq = C^. When we discuss leakage, we imply 
that the instantaneous state in C^ has elements in Z±i. In this case, the information within the state of the system 
cannot be described completely by the four operators above (Eq. ^). Since the total angular momentum with I = 
is a constant of the motion under the systcm-enviroment Hamiltonian, a state not completely represented by a linear 
combinations of (Eq. ^ is thus affected by decoherence. We will expore this DFS as implemented in liquid state 
NMR for both one and two logical qubits. 

The internal Hamiltonian (in the rotating frame) for two spins in liquid-state NMR already is exclusively in Zq and 
thus can be expressed by the operators in Eq. ^ it does not cause mixing of the subspaces Zi , 

H.nt = ^(<Ti-a2) + |Ji2ai.a^ (2) 

where Att'12 is the difference in chemical shift of the two spins and J12 the scalar coupling constant. The former 
coefficient scales the logical a^ operator, while the latter scales the a^ operator. Thus evolution under the inter- 
nal Hamiltonian alone generates a continuous rotation about an axis in the logical xz-plane making an angle of 
arctan(7rJi2/Act'i2) with the logical x-axis. As illustrated below, more general gates can be obtained via the interplay 
of the internal Hamiltonian and an external time-dependent RF (radiofrequency) field, with Hamiltonian 

H,,t = c.rf(i)e-^^*(*)Xe^^*(*) {X = ^E^O - (3) 

where (pit) is a time-dependent phase and (jJr^{t) is a time-dependent amplitude. The phase and amplitude are 
independently controllable. Note that -ffe^t cannot be expressed as a linear combination of the logical Pauli operators. 
In the presence of RF fields the evolution of a state inside the DFS under the total Hamiltonian necessarily causes 
the information to "leak" outside of the DFS, where it is no longer immune to collective dephasing. This combination 



of Hint and Hext{t) can generate any unitary in the C^ Hilbcrt space, guaranteeing universalit y. I n the absence of 
decoherence and assuming ideal controls, we can reach a unit fidelity for any desired gate |25l l26l|. Our interest is 
assessing the fidelity of control for a finite decoherence, for a finite bandwidth of our control parameters, and for 
Hamiltonians not respecting the symmetry of the logical subspace. 




Evolution time (units of tp, the n pulse time) 

FIG. 1: Shown above is the projection onto the logical subspace of a state initially inside the DFS, during application of 
an RF pulse for various ratios of -^. Defining the projection operator onto the logical subspace as Pl , we plot p{t) = 

Tr[{PLp{t)f]/Tr{p{tf), for t = ^ 2tp, where p{t) = e""*'rf'('"-+'"' Vf'e'"'rf*('"-+''') and uj,,tp = tt. The logical state completely 
returns to the subspace after application of a vr-pulse to both spins only when the spins have identical resonance frequencies 
{Au! = 0). If the ratio -^ is non-zero, as required for universality, the return to the logical subspace is imperfect (in particular, 
it is in general possible to go back to a state very close to the initial state in a time t > tp, but it is much more difficult 
to implement a n rotation). A lo gica l tt- pulse using a single period of RF modulation is not possible, a more complex RF 
modulation, like composite pulses |i27|. strongly-modulating pulses [la. [l7l| or optimal control theory p^. is required. In the 
above model, ^ — 500; the initial state of the system is a^ 



In the case of ideal control fields, an instantaneous 7r-pulse (tp — > 0) corresponds to a logical operation [llj |. since 
Px{7r) — e~*'^/^('^^+'^^) — — e~*'^/^('^^'^^\ which is equivalent to a tt pulse around a^. Figure Q] motivates the extent 
to which universality within the subsystem can be obtained in the finite tp regime. In this figure, we plot the purity 
of the projection of p{t) = 6"*'^'^*'^'^^+°'=: Vf'e*"'***^'^^"'''^:"-' on the logical subspace. In the limit of very high RF power 
(■^ -^ 0), the system undergoes a 7r-pulse in a time tp — ^^ and returns completely to the subspace after this time. 
It remains outside the subspace only for the duration of the pulse. For tOrf which are physically relevant (wi < 27rl00 
kHz and < Aw < 27r20 kHz), a single RF pulse does not result in a logical 7r-rotation due to off- resonance effects. 
Experimentally we are limited to finite tp and even our simple two logical qubit model system is sufficient to introduce 
several key challenges in implementing coherent control over logical qubits: (i) decoherence due to leakage outside the 
subspace during RF modulation periods, (ii) decoherence due to leakage outside the subspace after RF modulation, 
and (iii) loss of fidelity due to cumulative leakage with respect to the spectral density of the noise. 

Figure 121 shows an illustrative example of the integrated effects of a 7r-pulse applied to the two spins in such a DFS 
on the purity (Tr [p^] ) and correlation with the ideal final state (Tr [pwant/o]) as a function of the ratio of the relaxation 
rate I/T2 to the RF power ujrf- The initial states were chosen from the four logical Pauli operators, and we made the 
approximation that the internal Hamiltonian is zero during the application of these 7r-pulses. As would be expected, 
the desired result (negating the state in the case of po = ay,a^, or preserving it for po = l^i o-^ ) is rapidly degraded 
by the totally correlated decoherence during the 7r-pulse, unless the Rabi frequency is considerably faster than the 
relaxation rate. The increase in both the coherence and the correlation when the relaxation becomes fast compared 
to the rotation rate is due to a sort of "quantum Zeno" effect, so that the RF field itself is unable to rotate the state 
out of the DFS. In a complete analysis of the 2-spin case the effects shown in Fig. ^ must be combined to those in 
Fig. 13). 

Manipulating more than one logical qubit introduces further complexities to the control versus leakage problem. 
For the DFS considered, the extension to 2 logical qubits encoded into 4 physical qubits leads to the following basis 
states: 



|00)l^ 10101), 
|01)l^ 10110), 



|10)l^|1GG1), 
|11)l^|1G1G) 
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FIG. 2: Plots illustrating the loss of fidelity due to totally correlated decoherence during the application of a 7r-pulse about the 
a;-axis to the two spins of the DPS (see text). The dashed curves (red in the on-line version) are for the initial states po = 1 



or Po = a^ , while the lower curves (blue in the on-line version) are for po = o-y or po 



The left-hand plot shows the trace 



of p^ following the 7r-pulse as a function of the inverse product of the RP power tJrf and the relaxation time T2. The right-hand 
plot shows the correlation with the ideal final state, i.e. the trace of pwantP, following the 7r-pulse as a function of this same 
parameter. 



We define leakage any evolution that will cause the state to be not fully described by a linear combination of this 
four basis vectors. Since we focus on the challenges unique to controlling multiple logical qubits, we assume the 
internal Hamiltonian to be given by the Hamiltonians of each logical pair (as in cq. 0) and a coupling between two 
spins pertaining to two distinct logical qubits: 



H = Hli + Hl2 + ^71" J23 0^2 • i?3 



(5) 



The interaction term H2^ 
defined by the states 



0-2 



■ (73 couple the system initially in the subspace defined by the state Q to the subspace 



10011) and |1100), 



(6) 



-m/4//4 



for example e "/"'^s.-ilolOl) = iy=-|1100). If the noise is collective only over each pair of spins that encodes a logical 
qubit [23, the states © are not protected against it and will decohere. The internal Hamiltonian will be therefore 
responsible for leakage and the ultimate decay of the system. 

Notice that we would in general expect the noise to be collective over all the physical qubits, and not pairwise 
collective. In the case of NMR, this corresponds to a fluctuating external magnetic field, which is fully correlated. 
However, the differences in energies between qubits could be strong enough to effectively add a non-collective compo- 
nent to the noise. In particular, we can consider in NMR the case in which each pair is formed by spins of a different 
chemical species. In this case, the difference in gyromagnetic ratio makes the strength of the noise acting on each 
pair unequal, so that the noise is no longer collective. On the other hand, when the Zecman energy separation is 
considerable, the coupling between spins can be very well approximated by the diagonal part of i/j^g, i.e. cr^cr^, which 
does not cause leakage. 

When the noise generator is fully collective (as for homonuclear systems in NMR), the internal Hamiltonian still 
causes leakage, via the coupling to the states (©. Since these states belong to the zero eigenvalue subspace of the 
noise generator, they do not decohere. Information could still be lost at the measurement stage, since the states © 
arc not faithfully decoded to a physical states. A unitary operation is enough to correct for this type of leakage, and 
since decoherence is not an issue here, there are no concerns regarding the time scale over which the correction should 
be applied; however, amending for this unwanted evolution would in general mean the introduction of an external 
control, that, as seen, is a source of leakage leading to decoherence. 

For logical encodings other than the DFS considered, the natural Hamiltonian may drive the state out of the 
protected subspace even for single logical qubits; for example, the noiseless subsystem considered in Ref. |29l | will 
evolve out of the protected subspace whenever the chemical shifts or scalar couplings among its three constituent 
spins are not all equal. 

If we wish to do something more complicated than merely freeze the evolution of the system, e.g. to rotate the 
DFS qubits while simultaneously refocusing all the interqubit couplings, the complexity of the modulation sequence 
increases and the various causes of leakage will combine. In attempting to demonstrate a universal set of logic gates 



on a pair of two-spin DFS qubits by liquid-state NMR, leakage turned out to be an unavoidable problem for all 
practieal intents and purposes. Fortunately, it turns out that in many practical situations other means of inhibiting 
decoherence are also available, and can allow one to leave the protected subspace if need be in order to simplify the 
implementation of logic gates on encoded qubits. Dynamical decoupling is a particularly promising class of techniques 
for these purposes, which are applicable whenever the correlation time of the noise is long compared to the rate at 
which the system can be coherently modulated. The next section will analyze the principles involved in this approach, 
and show how they may be applied to some simple but realistic examples. 

III. STOCHASTIC LIOUVILLE THEORY AND CUMULANT AVERAGES 

The earliest example of a pulse sequence that could correct for random field fluctuations with long correlation times 
was given by Carr and Purcell pM l3lL l33 | . Today this would be regarded as dynamical decoupling or "bang-bang" 
control J33 and it has been applied beyond magnetic resonance, for example to the control of decoherence in spin- 
boson models p34| . In this section we outline a formalism, based on the well-known stochastic Liouville formalism 
|l5i I35I l3fil| and cumulant expansion, which allows us to analyze the effects of dynamical decoupling on decoherence. 
In the following section we apply this to the above two-spin DFS. 

Stochastic Liouville theory is based on a semiclassical model of decoherence, in which the Hamiltonian at any 
instant in time consists of a deterministic and a stochastic part. In the simplest case of NMR T2 relaxation, this 
typically takes the form 

Htotit) = Hdet{t) + Hst{t) ^ H™i + 77rf(t) + V,wfe(i)^/c , (7) 

^ — ^ k 

where Hint is the static internal Hamiltonian, Hrf{t) is the RF Hamiltonian, the ujk{t) describe the phase shifts due 
to stochastic, time-dependent fluctuating fields and Zk are the generators of each of these noise sources, i.e. operators 
which describe how these classical fields are coupled to the quantum system. In the two-spin DFS example considered 
previously, there is only one noise generator Z = {a\-\-a1)/2 with u){t) = 'yB{t), which describes collective fluctuations 
parallel to the applied static magnetic field. 

We now introduce a superoperator C(t) defined on Liouville (operator) space via 

£(i) = Htltit) ^ t - t ^ Htotit) = Cdetit) +Y,LUk{t)Zk (8) 

^ — ^k 

where Z^ — Z^ (^ TL ~ H (^ Zk- This superoperator is the generator of motion for density operator p, meaning 

p{t) = Up{Q) = rexp( -i I dt' £{t'yj p{0) (9) 

where T is the usual time ordering operator. Since what is actually observed in an experiment is the statistical 
average over the microscopic trajectories of the system {p{t)), we have to take the ensemble average superpropagator 
to obtain {p{t)) = (l/()p(0). The problem of calculating the average of the exponential of a stochastic operator has 
been solved by Kubo |37l | using the cumulant expansion. In terms of the so-called "cumulant averages" ( ■ ■ ■ )c (see 
Appendix A), the superpropagator is given by: 

(U) = cxp(-if dt'{C{t')),-\TJ dhf dt2{C{h)C{t2))c + ---] (10) 

Providing || L dt'C{t')\\ <C 1 for all i > we can safely neglect high order terms in the exponential's argument. 

Similar expressions are obtained in the formalism of average Hamiltonian theory (AHT) [Sg for the coherent 
(instead of stochastic) averaging of the system evolution under control Hamiltonians cyclic and periodic in time. Here 
we can obtain simplifications analogous to those encountered in AHT if we analyze the evolution in the interaction 
frame (called "toggling frame" in NMR) defined by the RF propagator Urf{t) [33. In this frame the noise operators 
acquire a further time-dependency (coherently imposed by the cyclic excitation) in addition to the stochastic time 
dependency of their coefficients ujk{t). The total Hamiltonian in the toggling frame is 

Htotit) = Hdetit)+Y,^i^kit)~Zkit). (11) 

where the toggling frame equivalent O of any given operator O is defined by 0(i) = U^^it) O C/rf(i), with : 

C/rf(i) = Texp(-iJ dt'H^it')\ (12) 



and Urf{tc) — 1 for cyclic controls, so that the toggling frame and laboratory frame coincide at the end of each cycle. 
This time-dependent change of basis in Liouville space induces a change of basis in the space of superoperators 
acting on Liouville space, as a result of which the noise super-generators Z^ also become time-dependent, i.e. 



Zkit) = Zk{t)(g,t - t(g>Zk{t) . 



(13) 



This facilitates the calculation of the average super-generator, and further allows the first-order effects of the RF 
fields upon the decoherence to be determined from the second-order terms in the cumulant expansion. In contrast, 
assuming as usual that the random variables Wfc(t) have a mean value of zero at all times, it would be necessary to 
analyze the third-order terms in order to obtain these results. 

Returning now to the problem of greatest interest here, in which there is only one noise generator which describes 
totally correlated decoherence as above and the corresponding random variable cj(<) is stationary and mean zero, the 
results given in Appendix A imply that the first two cumulants in the toggling frame are: 



JCiit) 



I [ dt'{Cdet{t')+Uj{t')Z{t')) = - [ dt'Cdetit') 

t Jo t Jo 



1 



(14) 



/C2(i) = 72 / dh 

t Jo Jo 



dt2([Cdet{tl),Cdet{t2)] + 2 G{t2 - h) Z{h) Z{t2) 



In the last line we have introduced the autocorrelation function G{At) = {uj{t + At)uj{t)) for the stationary random 
noise variable uj{t). 

IV. REFOCUSING NOISE WITH A CARR-PURCELL SEQUENCE 

We now use these results to analyze an implementation of a a^ rotation on a two-spin DFS qubit. We will show 
that this implementation is applicable when the correlation time of the noise Tc is long compared to the time required 
to apply a 7r-pulse to the spins . It consists of a (7r/2)-rotation of both spins in the DFS qubit about the y-axis, 
followed by a Carr-Purcell-style sequence consisting of an even number 2n of 7r-pulscs separated by equal time intervals 
T = t/2n, and finally the inverse (7r/2)-rotation, i.e. 



H. 



(15) 



This transforms the weak (Tz<Tz coupling between the two spins of the DFS qubit into a-xf^x, which projects to the a^ 
operator within the DFS (Eq. Q). Setting r = (f)/{2n-Kj) thus yields a rotation by an angle 4> around the logical 
X-axis. Even though the state of the two spins is outside the DFS throughout the time 2nT, the sequence of 7r-pulses 
is able to refocus the effects of the noise provided t <^ Tc- 

Assuming instantaneous 7r-pulses, this follows from AHT since during any cycle (0, 2t) the internal Hamiltonian in 
the togghng frame Hint alternates between -|-Atj((T^ — a^.) + (7r/2) Jct^ct^ (in the interval (0, r)) and — Aijj(crj|, — cr^) -|- 
(7r/2) Jcr^cr^ (in the interval (t, 2r)), so that the zcroth-order average Hamiltonian is just H'^'^^ = {-k /A)al.a'^ = (7r/4)cr^. 
This is in fact also the average Hamiltonian to all orders, since the toggling frame Hamiltonian commutes at all time, 
and the first cumulant is just the corresponding superoperator K,i = K.i = {2t)^^{H* ^l — t®!!). 

Again because the toggling frame Hamiltonians commute, the deterministic part of the Liouvillian Cdet{t) does 
not contribute to /C2 = /C2 at the end of each cycle, nor at the end of the entire sequence. The second cumulant is 
therefore determined by the stochastic part alone: 



IC2 



(2?it) 



2nr f^ti 

dti / dt2G{t2-h)Z{ti)Z{t2) 

Jo 



(16) 



Because each 7r-pulse simply changes the sign of Z{i) from the preceding interval, it follows that Z{i) = +Zx if t is 
in an even interval (2fcT, (2fc -I- 1)t) and Z — —Zx if t is in an odd interval [{2k — l)r, 2fcr), where k is an integer 
< k < n and Z^ is the noise super-generator rotated along the x-axis. In addition, since the random variable w(t) 
is stationary, the double integral over any two intervals i and i' will be equal to the double integral over any other 
pair J and j' providing that \i — i'\ = \j — j'\. These observations allow the overall double integral in Eq. (|16|l to be 
expressed as 



JC2 



2Z2 

(2nT)2 



2nA+ ^ {2n-m)Br, 



ZlC 



(17) 



where 



and 



A = 



/ dh / dt2G{ti-t2) 

Jo Jo 



B,„ = (-1)" / dii / dhG{ti-t2) 

JtriT Jo 



for m = 1, . . . ,2n — 1. 

In the case of Gaussian noise with autocorrelation function G(t) = Q'^e^*'/^", one finds that 



A 



{nr^y (e-^/^" + t/t, - 1) 



and 



On evaluating the geometric series in Eq. H17|l . one obtains the closed form 



c 



2rj2 



(2nT)2 



2n( t/tc + e 



— r/r 



1 



1 - e-^/^ 
1 + e-^/^ 



1 -2n(l + e 



1)^ 



— t/Tc 



^-2nr/ 



'■•) 



(18) 
(19) 

(20) 
(21) 

(22) 



which is easily shown to go to zero as t/tc — > [43| ■ We can quantify the protection afforded by the CP-sequence by 
taking the entanglement fidelity |l3l |40j of the superoperator with the ideal propagator for the sequence as a measure 
of its efficacy, F = Tr [f/j^^S*] . Since the unitary part of the evolution commutes with the noise and gives the ideal 
propagator, the fidelity is just the trace of the superoperator, which for a single two-spin DFS qubit is 



^(C) 



Tr 



-Z2c(2«t)V2 



|(3 + 4e-2C»'^'+e-'«"'"') 



(23) 



The fidelity for cycles of CP-sequences of length 4 and 16 are plotted in Fig. 3. As expected, it shows an improvement 
for an higher number of intervals and shorter time spacings with respect to the correlation time. 
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FIG. 3: Gate fidelity as a function of the correlation time for 4 and 16 cycles of the Carr-Purcell (CP) and Time-Suspension 
sequences (TS). The noise strength il was fixed at 1 Hz., while the duration of the entire sequence was fixed at ttot = 4 sec 
(where ttot ~ 2nr for the CP sequence and ttot ~ 4nr for the TS sequence). The increase in fidelity at very short correlation 
times is due to the phase fluctuations becoming so fast that they produce essentially no effect at the given noise strength f2. 

It is interesting to also consider a simple sequence that completely refocuses the internal Hamiltonian, namely 

i---[-]i---[-]i---[-]i---[-]iy^ (24) 



where the superscripts on the pulse angles now refer to the spin afFeeted by the pulses and n is an integer. This will 
be referred to in the following as the Time-Suspension (TS) sequence. The average Hamiltonian is now zero, while the 
noise operator in the toggling frame is Zi = ±(cr] +o-^)/2 in the intervals 1 and 3 respectively and Z2 = ±(0-] — (t^)/2 
in the other two intervals. If we sandwich the TS-sequence between a pair of (7r/2)-pulses as we did for the CP, 
and again assume a stationary and Markovian Gaussian distribution of totally correlated noise, we find it is more 
effective at protecting the system from decoherence even when the number of 7r-pulses on each spin and the cycle 
time is the same, since the effective modulation rate is then faster (there is a pulse every t/2). Indeed the relaxation 



superoperator for the TS-sequence is IC2 
{k = 1,2) and: 



Ci{2f + Z2) + C2 ^1^2 (see Appendix B), where Zk 



1-1 



Ci 



16n-^T^ 



l-e-/- 



^—Anr/Tc 



^gir/T^ _ ^ 



+2n^^ + n (e-2^/^- - l) (2 - fi-^/^-) 



1) - 1) e-3^/^<= 



(25) 



C2 = 



1 



oT/Tc 



-4T/Te 



167i2t2 1 + e2^/ 

The fidelity is therefore: 



-Anr/- 



"-" [ne^^'^" -n+l\+ ne*^/^'= -{n+l) 



Tr [exp {-Ci{Zl + Z|)(4nT)V2 - C22i22(4nT)V2)] 
1 e-Ci(4"T)^ (^ cosh(-Ci(4nT)2) -f cosh(-C2(4nT)V2)) 



(26) 



(27) 



V. SIMULATION OF A SELECTIVE DPS QUBIT GATE 



The analytical expressions found above for the attenuation due to totally correlated noise with a stationary Gaussian 
Markov distribution apply only to the special case of ideal pulses (instantaneous in time), but similar behavior is 
expected under more realistic assumptions on the control fields. In particular, to act selectively only on some of the 
spins we would have to use the technique of SMP |l6l llTj , thereby inducing a much more complex dynamics on the 
system for which closed form solutions are not available, but which can be studied via numerical simulations. 

Ci,C2 I " I " I ^ I ^ II 

Ti: 71 Ti: 71 Ti-n 

cLj i_^ 



71; 



n 



n 



n 



FIG. 4: Selective rotation about the logical a;-axis of a two-spin DFS qubit, while the evolution of a second DPS qubit under 
the internal Hamiltonian of the system is refocused. 



We have studied the accuracy with which a rotation about the logical x-axis can be performed by numerical 
simulations. These simulations included the internal Hamiltonian, the external control Hamiltonian and totally 
correlated noise uj{t) with a stationary, Markovian Gaussian distribution. The evolution was discretized into equal 
time steps, for each of which we calculated the propagator C/(ifc) = exp(—i(Hi nt + Hrf{tk) + ^{tk)Z)St). The noise 
strength uj{tk) is extracted from a multivariate gaussian probability distribution |4J|j with a covariance matrix Cj fc = 
Sl2e~IJ~'=l'5*/'rc^ where j and k are integers indicating the time intervals. We then take the average of the superoperators 
Si = Ui® Ui obtained over a sequence of evolutions differing only by the random number seed. 

We have performed one set of simulations using a fictitious two-spin molecule (chemical shift difference: Alo = 
600Hz, scalar coupling J = 50Hz), and another using the internal Hamiltonian of ^"^C-labeled crotonic acid, a 
molecule containing four carbon spins [41|. Both sets of simulations were performed with instantaneous ideal pulses, 
and again with the strongly- modulating pulses used in actual NMR experiments. SMP are time-depend RF fields 
designed by a numerical search, and perform precise rotations of one or more spins while refocusing the evolution of 
all other spins in a molecule [la. Il7| . 



In the case of the two-spin molecule, since selective pulses are not required, we compare the results of SMP pulses 
with the dynamics under short, collective pulses (called "hard pulses", Tr-pulse time tp = 2/xs). SMP appear to 
perform better even if they require longer times. In the crotonic acid simulations, the sequence was designed not only 
to implement a selective 7r/2-rotation about the logical x-axis on the two spins in one DFS qubit, but to also refocus 
the evolution of the other two spins under the molecule's internal spin Hamiltonian (see Fig. 4). 

The fidelities of these simulations are plotted as a function of correlation time in Fig. 5. Compared to simulations 
with ideal pulses, we observe a drop in the fidelity due to the finite duration of each pulse. This drop is only in 
part accounted for by the increase in time in the cycle length. Nevertheless, the effectiveness of the CP-sequence in 
preventing decohcrcncc during the unavoidable excursions from the DFS is evident. 
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FIG. 5: Fidelity for ideal and real pulses, a) Two-spin CP sequence implementing a 7r/2 rotation about the logical a^ (A 
fictitious spin system with Auj = 600Hz and J=50Hz was used in the simulation) b) tt/2 rotation about the logical a^ for 
Crotonic acid (see Ref. 1411 for chemical shifts and J-coupling data). 



VI. CONCLUSIONS 

In this paper we have considered the difficulties of operating on quantum information stored in encoded qubits with- 
out losing the protection from decohcrcncc offered by the encoding. Although we have focused upon the Hamiltonians 
and control fields operative in NMR for concreteness, similar difficulties will be encountered in other functional realiza- 
tions of quantum information processing, including squids, ion traps and quantum optics. The most significant result 
is a demonstration that in many realizations, including NMR, the implementation of a universal set of quantum gates 
may be considerably simplified by briefly leaving the DFS while using dynamical decoupling to inhibit decoherence 
during these excursions. This approach depends on the ability to operate on the system on time scales short compared 
to the correlation time of the noise. In evaluating various possible realizations of quantum information processing, it 
is important to characterize not only the decoherence rate, but also the spectral density of the underlying noise, to 
verify that the gate speed is sufficient to allow the noise to be refocused. In addition to its role as a facile testbed for 
quantum information processing, NMR spectroscopy provides widely applicable tools such as Average Hamiltonian 
theory by which one can calculate the efficacy of control sequences for rcfocusing the noise, and devise new ones for 
specific noise generators. 

Acknowledgments: This work has been supported by the US Army Research Office under grant number DAAD19-01- 
1-0519, and by the Quantum Technologies Group of the Cambridge-MIT Institute. 



VII. APPENDIX: THE CUMULANT EXPANSION 



We calculate the ensemble average of a time-ordered exponential in terms of the cumulant expansion. First, expand 
the time-ordered average exponential S = (Texp(— i J„ dt'H{t'))) via the Dyson series [i^ : 
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(28) 
+ ^^/o dti--- /o dtn{H(h) ■ ■ ■ HiU,)) + ■■■ 

The term (iJ(ii) • • • H(tn)) is called the n-th moment of the distribution. We want now to express this same propagator 
in terms of the cumulant function K{t), defined by: 

S = e^(*) (29) 

The cumulant function itself can most generally be expressed as a power series in time: 

K{t) = V ^^A'„ = -itK, + ^-^K, + ■■■ (30) 

^^^ n\ 2! 

Tl = l 

Expanding now the exponential (|29|1 using the expression in equation l|3()() we have: 

S = l + K{t) + ^{K{i)f + --- 

where in the second line we have separated terms of the same order in time. By equating terms of the same order in 
Eq. (|31|) and H28|) we obtain the cumulants Kn in terms of the moments of order at most n. For example: 

K\^ \j\t'{H{t')) 

(32) 



K2= ^T J dh J dt2 {H{h)H{t2)) - Kl 



The propagator can therefore be expressed in terms of the cumulant averages: 

{H{t'))^ = {Hit')) 
{H{t,)H{t2))^ = T{H{t,)H{t2)) - {H{t,)){H{t2)) 

The propagator can therefore be written as: 

S = exp(-i f dt'{H{t'))^ - I dti f dt2{H{h)H{t2))^ + ---) (34) 

Note that if i? is a deterministic function of time, the ensemble averages can be dropped and (iJ(t)) = /„ dt' H{t') 
becomes the time-average Hamiltonian, which is the first term in the Magnus expansion. The second term in the 
cumulant expansion, on the other hand, becomes 



T I dti I ' dt2H{ti)H{t2) -{ I dt' H{t' 
Jo Jo \Jo 

^ 2 [ dti r dt2H{ti)H{t2)- [ dti [ dt2H{ti)H{t2) 
Jo Jo Jq Jo 

[ dti r dt2H{ti)H{t2)- [ dti [ dt2H{ti)H{t2) 
Jo Jo Jo Ju 

dti / dt2 [-ff(tl), H{t2)] 

Jo 



where [•,•] denotes the commutator and we have used the fact that the time-ordering operator T symmetrizes its 
argument with respect to permutation of the time points. This is the second term in the Magnus expansion for the 
"average" (effective) Hamiltonian, and proceeding in this fashion one can in principle derive average Hamiltonian 
theory from the Dyson and cumulant expansions. 
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VIII. APPENDIX: CUMULANTS OF THE TIME SUSPENSION SEQUENCE 

We derive equations 1251 and 1261 for the fidelity attenuation of a Time Suspension sequence. Consider first the basic 
cycle of the sequence, composed of four time intervals. The average internal Hamiltonian is now zero, while the noise 
operator in the toggling frame is ±Zi = i{crl + a1) in the intervals 1 and 3 respectively and ±^2 — ilc^ — a^) in 
the other two intervals. 
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FIG. 6: Domain of integration and toggling frame noise operator for the TS sequence 



By inspection of the domain of integration, the second order cumulant for the first cycle, (corresponding to the first 
triangle in figure O is given by: 



K. 



(t) 



^(2A - Z2B3,i) + (1% (2A - ^4,2 



(4t) 
(4r)2 



(2A-Be~2^) 



^^B(e- 



■^5pp-(S2,l - ^4,1 - -83^2 
-3^ 



^4,3) 



(4t 







(36) 



where A and B where defined in equations H2U|) and (|21|) and B^^k = Be^'-'^^''-''^/'^''. If the sequence is repeated 
n-times, we can divide the domain of integration in 4 by 4 time step squares and triangles (where the triangles are 
equivalent to the first cycle and the first of the possible squares is the shaded area in figure El . Each triangle will give 
the same contribution calculated above for the first cyclc(if there are n cycles, we will have n of them). The second 
cumulant from the first 4 by 4 square, corresponding to ti in the second cycle and t2 in the first cycle, is: 



IC. 



is) 



7-2 I 72 

^1 + ^2 o/r,„-4^ 



(4r)2 



-S(2e" 



-2^ 



-6—\ , ^1^2 n^ -3 — 



(4r) 



+ e 



'5^ 



e''^) 



(37) 



-(^) 



To calculate the contributions from ti in cycle k and ^2 in cycle h, it is enough to multiply this cumulant IC2 by 
e~ ~^'~ '. In general, for n cycles we obtain: 



fc-i 



IC2 



- (ri/C- 



(t) 
2 



k=2h=l 



-i^ik-h)^(s) 



/C^^0 = Cl(^l+^2)+C2^1^2 



(38) 



with: 



16n2r2 



(J^^^)' (e-4-/r. („g4r/r. _ ^ + l) - l) g-^^/^^ 



+2n:^ + n (e^^^/^- " l) (2 - e"^/^<= 



(39) 



and 



16n^T^ 



fl - e^/^= 



-4r/r, 



1 + e2^/^= 



^— 4nr/T 



ne^^/^" - n 



a'tr/r^ _ 



(n + l) 



(40) 
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